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Abstract. I solve a Relativistic Navier Stokes Model assuming boost invariance and 
rotational symmetry. I compare the resulting numerical solutions for two limiting 
models of the shear viscosity. In the first model the shear viscosity is made proportional 
to the temperature. Thus, r\ oc T/ao where <7q is some fixed cross section (perhaps 
do ~ Aq CD ) . This viscosity model is typical of the classical Boltzmann simulations 
of Gyulassy and Molnar. In the second model the shear viscosity is made proportional 
to T 3 . This model is typical of high temperature QCD. When the initial mean free 
path of the T 3 model is four times larger than the T/ao model, the two models of 
viscosity produce the same radial flow. This result can be understood with simple 
scaling arguments. Thus, the large transport opacity needed in classical Boltzmann 
simulations is in part an artifact of the fixed scale ao in these models. 

1. Motivation: The observation of elliptic flow is one of the most striking results 
the heavy ion program X\. In mid-peripheral collisions the azimuthal anisotropy of 
the produced particles ^(pt) = ( COS (20)) PT rises linearly as a function of transverse 
momentum up to pt ~ 1.5 GeV and then flattens and maintains a constant value of 
approximately 15%. Ideal hydrodynamics provides an economical description of the 
observed elliptic flow for transverse momentum below 1.5 GeV [2]. However, ideal 
hydrodynamics assumes that the transport mean free path £ m f p path is so small that 
viscous terms can be neglected. Simple estimates indicate that that in heavy ion 
collisions this assumption is marginal. It is hoped that £ m f p /L is of order ~ 1/5 and 
that some semblance of an equilibrium QGP plasma will be formed. For £ m f p /L ~ 1/5 
viscous hydrodynamics should provide a semi-quantitative guide to the evolution of the 
system. Viscous hydrodynamics should then explain the systematics of the elliptic flow 
measurements. 

Elliptic flow was studied in the fully non-equilibrium framework of classical kinetic 
theory by Gyulassy and Molnar These authors varied the cross section between 
classical particles and computed the resulting elliptic flow. The results are surprising. 
First, the full kinetic theory reproduces the full shape of the V2{pt) curve. In particular 
the kinetic theory describes the linear rise with transverse momentum (as expected from 
hydrodynamics) and the subsequent flattening of ^(pr) • The first viscous correction 
to the thermal distribution function clarifies this transition from the hydrodynamic to 
the kinetic regime j3j. In spite of this success, these classical Boltzmann simulations 
suggest that hydrodynamics is not responsible for the observed t>2- The cross sections 
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needed to reproduce the observed elliptic flow are « 10 mb. With such cross sections 
the initial momentum degradation length is approximately four times smaller than the 
thermal wavelength. This seems impossible. 

However, several differences between between QCD and these classical cascade 
models should be noted. First, the classical cascade preserves particle number. Second, 
the classical simulation have a fixed scale ctq. This is not so unreasonable - perhaps o~q 
is A.q CD . To understand the implications of this fixed scale let us estimate the temporal 
dependence of the viscosity in different situations. 

In the perturbative quark gluon plasma the shear viscosity is proportional T 3 as 
given by dimension. Thus the ratio of the shear viscosity to the entropy rj/s is a 
constant up to logarithms. In a classical massless gas with constant cross sections the 
shear viscosity is given by rj = 1.264 T/<t Thus the temperature dependence of the shear 
viscosity in the classical model j3] differs from high temperature QCD. Simple scaling 
arguments indicate that a model with a constant cross section is unlikely to thermalize. 
For a Bjorken expansion, the condition for hydrodynamics to be valid is: Y s /t <C 1 
where where T s is the sound attenuation length p]/{e + p) . 

Next, I estimate how T s /r evolves as a function of time during a Bjorken expansion 
for the two limiting models of viscosity discussed above. For a Bjorken expansion 
T oc 1/r 1 / 3 and n oc 1/r. Then, when the viscosity is proportional to T 3 we find 
that the system comes closer to equilibrium as a function of time 
T s 1 1 

— ~ — = ~ -ttr (ID Bjorken with ri oc T 3 ). 

T Ti T 2 ' 6 

In contrast for a classical gas with conserved particle number e + p = 4nT and constant 
cross section the degree of thermalization remains constant as a function of time 

r s 1 

— ~ ~ Const (ID Bjorken rj oc T/o~q) . 

t rna 

Similar instructive arguments may be given when the system expands in three 
dimensions. The conclusion remains the same. When the viscosity is proportional 
to T 3 the system is much more likely to thermalize than when other scales such as <Jq 
and Aqcd enter the problem. 

2. Viscous Solutions. It has been understood for some time that the relativistic 
Navier stokes equations can not be solved directly. However, this problem is cosmetic 
rather than fundamental in nature 0. A myriad of different hydrodynamic models [HIE] 
can be solved numerically with considerable effort jH]. These models all give the same 
solution up to corrections which are proportional to (T s /t) 2 . The stress energy tensor 
is always close to its canonical form T ,J ~ r/ [d % v^ + &*v l — |5 U div 1 ^ . These model 
equations have no greater validity than the Navier Stokes equation. The hydrodynamic 
model used here was inspired by a model with appealing mathematical structure due to 
Ottinger [H]. 

First, I solve a Navier Stokes Model (the Ottinger Model 6\) assuming boost 
invariance and radial symmetry. The initial conditions are taken from P. Kolb's inviscid 
hydrodynamic calculations which reproduces the observed spectra and multiplicity. 
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Figure 1. A comparison of the viscous and inviscid (Euler) solutions for a 
Bjorken expansion with radial symmetry. The initial conditions are from P. Kolb's 
hydrodynamic model The equation of state is p — \e. The viscosity is proportional 
to the entropy r\js =\. (a) The energy density multiplied by r for various times (b) 
The fluid velocity for various times. 



The viscosity is made proportional to the entropy rj/s = 1/5. The viscous solution is 
compared with the inviscid (Euler) solution shown in Fig.^ Examine the energy density 
in Fig. I^a). First the viscous solution does less longitudinal work since the longitudinal 
pressure is reduced by the longitudinal expansion (see e.g. |3j). Consequently, the energy 
density initially decreases more slowly for the viscous case. However, the transverse 
pressure is increased by longitudinal expansion. This causes the transverse flow to rise 
more rapidly in the viscous case as seen in Fig.QJb). This larger transverse flow velocity 
subsequently causes the energy density to fall more rapidly in the viscous case. By a 
time of ~ 6fm the viscous and inviscid solutions are similar. In summary, viscous 
corrections do not integrate to yield an order one change to the inviscid flow. 

Next, I compare two simple models for the viscosity. The first model for the viscosity 
is taken from a classical ideal gas with a constant cross sections. In this case r\ = 1.264— 
with cr = 10 mb. This model of the viscosity has been studied within the domain of 
kinetic theory by Gyulassy and Molnar 

The second model for the viscosity is referred to as the Minimal Model below. In 
this model we take 

| 1-264 ^ for e < e c 
1 | s for e > e c 

where e c — I GeV/fm 3 and o"o = 10 mb. This model of the shear viscosity has rj cx T 3 
for high temperatures but has a fixed scale do (i.e. Aq^ d ) at low temperatures. The 
shear viscosity of the Minimal Model is always larger than the fixed cross section model. 
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Figure 2. A comparison of two models of the shear viscosity. In the first model the 
shear viscosity r\ = 1.264T/<To. The minimal model of viscosity is described in the text 
and has ij oc T 3 when the energy density exceeds 1.0GeV/fm 3 . In both case p = ^e. 
The thin lines are lines of constant transverse rapidity with values 0.1, 0.2, 0.3 ... The 
thick line labeled Freezeout is where the viscous "correction" is half of the pressure. 
(dt) denotes the entropy weighted mean transverse rapidity along the freezeout curve. 



The viscous hydro dynamic solutions to these models are illustrated in Fig. El The 
two models of viscosity give approximately the same solution although the viscosity in 
the Minimal Model is initially four times larger than the fixed cross section model. 
Eventually, the viscosity becomes large and the system freezes out. The curve 
labeled "freezeout" indicates when the viscous correction becomes equal to half of the 
hydrodynamic pressure. The mean transverse relativistic velocity (ut) is calculated 
along the "freezeout" curve using the Cooper- Frye formula and weighting each surface 
element by the entropy. The two models of the shear viscosity give approximately the 
same transverse flow as can be seen be comparing [yj) in each case. Thus the large 
transport opacity needed in classical Boltzmann simulations is in part an artifact of the 
fixed scale o~o introduced into the problem. 
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